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We study a coupling flow of pure QCD gauge system by using the Monte Carlo Renormalization Group method. 
A rough location of the renormalized trajectory in two coupling space is obtained. Also we compare 4 different ac- 
tions; (a)standard Wilson, (b)Symanzik's, (c)Iwasaki's and (d)QCDTARO's. The rotational symmetry is restored 
better as an action gets close to the renormalized trajectory. 



1. Introduction 

In order to achieve the continuum physics in 
lattice gauge theory, a lattice spacing a should be 
small enough. However the scaling violation of 
Wilson action at presently accessible a is clearly 
seen fl]-||] . Since the cost of a Monte Carlo Sim- 
ulation of QCD can estimated as oc a~ 6 it 
seems difficult to adopt smaller a without a sur- 
prising computer power- up or improvement of the 
algorithm. 

Recently "improved" actions have been exten- 
sively applied to resolve the scaling problems j| . 
Although in these studies, several types of actions 
are proposed and used. What is the most opti- 
mal choice or what is a "perfect" action is still an 
open question. 

In this paper, we analyze the coupling flow 
of pure SU(3) gauge theory by the Monte Carlo 
Renormalization Group (MCRG) method. And 
we compare several actions, which are widely used 
in actual simulations, by investigating the rota- 



tional symmetry. 

2. Improved Action 

In this paper we consider 2-coupling space in 
which the lxl simple plaquette and 1x2 rect- 
angle loop are included. The action S is, 

S = /3n^ReTr(l-~!7 p/og ) (1) 

plaq 

+ /3i2^ReTr(l-if/ rect ), 

rect 



where 




Although there are several arguments about 
the necessity of other loops as possible candidates 
of an action |||| , actions with two couplings are 
widely used as improved ones in actual lattice 
QCD simulations Q. Here, we restrict ourselves 
on a 2 coupling space and our purpose of this 
paper is to propose an optimal choice of action in 



2 coupling space. Inclusion of other couplings is 
now in progress and will be presented elsewhere. 

For actions with 2 couplings, we consider 4 can- 
didates which are widely used; (a) Wilson || , (b) 
Symanzik §, (c) Iwasaki §, (d) QCDTARO §]. 

/3s (or the ratio 712 = finl in these actions 
are listed in Table |[ 

Table 1 

couplings in each action 
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(a) Wilson 









(b) Symanzik 


-0.05 






(c)Iwasaki 


-0.091 






(d)QCDTARO 




6.1564(53) 


-0.6241(23) 






7.986(13) 


-0.9169(41) 



3. Coupling Flow by MCRG 

First we follow the coupling flow by MCRG. As 
MCRG method, here we adopt Swendsen's factor 
2 blocking scheme (9). For details of this blocking, 
see §,!. 

3.1. Schwinger-Dyson method for determi- 
nation of effective action 

From the original configurations {U} which is 
generated by an action S with (/3n, P12) in eq.(^) 
we obtain the blocked ones {U 1 } after the above 
blocking procedure. When the blocked configura- 
tions {{/'} can be considered as generated config- 
urations by an action S' with ((3'u, f3' 12 ) m ec f(H)! 
an action S' is an effective action of the blocked 
configurations, and 



(2) 



can be regarded as the coupling flow associated 
with this blocking. 

For the determination of an effective action, 
here we use a Schwinger-Dyson method [|Io|] . 

This method is simply based on the following 
identity. For a link Ui, we consider a quantity; 



F = IT / dUilmTr{X a UiGf) 



s 2ReTr(£/ i G ! ) 



(3) 



where I means link except for I. G\ is a sum 
of staples Gf for link / as G; = J2 a § G T- For 



infinitesimal transformation; U\ — » (1 + ie\ a )Ui, 
invariance of F leads to 

^ = YljdUjlReTrdX^UtGr) (4) 
+ ImTr(A Q [/;GnimTr(A Q [/;G0] 

x e 2ReTr(C/iG,) = q 

Using a well-known formula of A-matrix and 
summation over a and integral over link Ui , we 
have 

^Re < TrtUtG?) > = (5) 

Y,^{-Re<Tr(U l G?U l G])> 

7 

+Re<Tr(Gr(G/) t ) > 

--Re < lmTr(UiGf)ImTr(UiG]) >} . 

We then apply this equation to the blocked con- 
figurations. 

3.2. Lattice Simulation and Results 

We use 8 4 lattice size and about 2000 config- 
urations separated by every 10 sweeps are used. 
For the starting point (/3n,/3i2), we try (a)-(d) 
coupling sets in Table [l] and other points are also 
tested to search the renormalized trajectory(RT) 
in detail. All simulations have been done on 
CRAY J90 at Information Processing Center, Hi- 
roshima University and on VPP500 at KEK (Na- 
tional Laboratory for High Energy Physics). 

In Fig.l, the result of coupling flows in 
(Pll, P12) plane is shown. Arrow denotes the mea- 
sured coupling flow in eq.(|^). 

Ending points seem lie on the universal line 
which can be thought as the RT of this blocking 
scheme. For comparison, we prepare 2 sets of (a)- 
(d) points at almost the same scale; a 2 a ~ 1.51(6) 
and a 2 a ~ 0.82(6). Among (a)-(d) at almost the 
same scale, (d) is the closest to this RT and (c) 
is the next. 

4. Rotational Invariance 

In order to see the improvement of the action, 
here we check the rotational invariance by mea- 
suring the heavy quark potential. In this paper, 
we compare 4 actions in Table [j]at a 2 a ~ 1.51(6). 



In FigJ^, the result is shown. Since the lattice 
spacing at this scale is rather large, rotational 
symmetry can not be expected to be restored on 
the standard Wilson action, which is clearly seen 
in FigJ|(a). For Symanzik action in Fig.||(b), 
the situation is similar to Wilson case. On the 
other hand for Iwasaki and QCDTARO cases in 
Figs.||(c) and (d), the rotational invariance can 
be seen. It is also noted that errors become small 
as the action gets close to RT. 

In order to compare quantitatively, we fit the 
measured potential V(R{) (Ri is the measured 
point) with the linear plus coulomb term as 
f(R) = a + bR + c/R and evaluate the effective- 
ness of this fit by 



x 2 = J2( v (Ri)-f(Ri)f 



(6) 
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QCDTARO 
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From these results, rotational symmetry can be 
restored better as the action gets close to RT. 
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Figure 2. Heavy quark potential at about a 2 a ~ 
1.51. The circles are measured values and dots 
are fitting values, for (a) Wilson, (b) Symanzik, 
(c) Iwasaki and (d) QCDTARO. 
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Figure 1. Coupling Flow in ((3n,(3i2) Space 



